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In this work it is shown, that for short 3ns neutrino pulses reported by OPERA, a relativistic 
shape deforming effect of the neutrino distribution function due to spontaneous emission, produces 
an earher arrival of 65.8ras in agreement with the reported 62.1ns± 3.7ns, with a RMS of 16.4ns 
explaining the apparent superluminal effect. It is also shown, that early arrival of long 10500ns 
neutrinos pulse to Gran Sasso, by 57.8ns with respect to the speed of light, could be explained by 
a shape deforming effect due to a combination of stimulated absorption and spontaneous emission, 
while traveling by the decay tunnel that acts as a LASER tube. 

I. INTRODUCTION 

Recently, OPERA collaboration has reported 20 superluminal neutrinos coming from 3ns LHC extractions, in 
addition to the 16000 neutrinos, collected for 3 years coming from 10500ns LHC extractions [1] that apparently are 
faster than light for w 60ns while traveling 730 km from CERN to Gran Sasso National Laboratory. What is most 
amazing is that the velocity obtained is even faster than the known speed of light in the vacuum, within margins of 
error accurately calculated. This result is in direct contradiction with Einstein's Relativity Theory that has been the 
angular corner of modern physics for more than a century. The speed of light in vacuum is considered to be the limit 
velocity for all matter and information travel, as a consequence of the causality principle. 

OPERA faster than light neutrinos is in direct contradiction with astronomical observations as the SuperNova 
1987A in 1987 three independent neutrino observatories: Kamiokande II in Japan, 1MB in USA, and Baksan 
in the former USSR detected 11, 8 and 5 neutrino events, in a burst lasting less than 13s almost at the same UTC 
time. Approximately three hours later the light from the Nova was observed. This does not indicate that neutrinos 
arrived faster than light, the accepted explanation is that neutrinos arise from the collapse of the star core, but the 
burst of light occurs only when the shock waves reach the star surface. If OPERA result is right, the neutrinos must 
precede the light by 60ns for each 730A;m, as the Nova 1987A is 168.000 light year away from the earth, neutrinos 
must arrived earth 1500 days BEFORE the explosion was observed, an not few hours earlier. 

Even recently, it has been argued that superluminal neutrino will decay very fast by Cherenkov analogue effect due to 
neutral current inter actions [5|]- ICARUS collaboration, another team at Gran Sasso, had reported that superluminal 
CNGS neutrinos do not decay as theoretically expected questioning OPERA result. To my point of view this 
objection to OPERA results is not valid because is based in a relativistic model, that will be not longer valid if 
superluminal neutrinos do exist. It is the same as to said that Bohr atom is forbidden because of Larmor formula will 
make the electrons to radiate. 

In a recent work [?], it has been proposed an explanation for early arrival by 57.8ns of the neutrinos with respect to 
the speed of light coming from "long" 10500ns LHC extractions. This could be explained by a shape deforming effect 
of the neutrino distribution function, with respect to the proton distribution function (PDF) due to the stimulated 
absorption. But in this work it was also probed that stimulated absorption can not explain the apparent earlier 
arrival of 20 neutrinos reported for "short" 3ns LHC extractions. The objective of this work is to probe that the 
shape deforming effect due to spontaneous emission produces a backward shift in time of 65.8ns theoretically calculated 
on the basis of special relativity LASER equations. That is in agreement with the reported result. 

II. THE EXPERIMENT 

First of all it is necessary to outline the experiment. The proton beam is produced with the CERN Super Proton 
Synchrotron (SPS), these protons are ejected with a kicker magnet, in two extractions each lasting 10.5yLts and 
separated by 5077is towards a graphite target where billions of mesons are produced. These mesons are focused into a 
1km vacuum tunnel where the mesons decays into muons and neutrinos. Then neutrinos continue traveling through 
the inside of the earth by liQkm until they arrive to Gran Sasso Laboratory 2.4ms later. Neutrinos are detected by 
OPERA in two separated groups: the first with mean neutrino energy of 13,9GeF and the second with A2,9GeV, 
corresponding to each one of the two proton extractions. 

Proton extractions are similar to step functions with several peak or oscillations superposed. These peaks cor- 
responds to the proton synchrotron radio frequencies of the SPS and the kicker magnet. The form of the time 
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distribution of proton function is accurately measured by a fast Beam Current Transformer (BCT) at center of the 
graphite target. The key of the experiment is that each maximum of the neutrinos detection must correspond to a 
maximum of the proton intensity. 

To obtain the time travehng of the neutrinos, each of 16000 neutrino events detected were tag in time, and correlated 
with its corresponding proton time distribution function for each extraction with high accuracy. To do that, a 
probability density function (PDF) is constructed, summing up all the proton time distribution function, for which 
neutrino interactions were observed at the detector. Then this function is shifted in time (t — > i + TOFc) by the 
estimated time of flight {TOFc = 730.085 fcm ^ ^j^g speed of light c. The peaks of these PDF function must correspond 
in time with the peaks in the neutrino detection if they fly exactly at the speed of light. The measured neutrino time 
distribution, detected at OPERA must have a delay time {t ^ t + TOF^), corresponding to the time of flight of the 
neutrinos {TOF^ = 730.o^5fcm ^^^^ velocity v. As both time functions, the theoretical PDF function and the 

detected neutrinos distribution had been shifted by TOFc and TOF^, respectively; the maximum likehood analysis 
must gives the best fitting for a time lapse [l|: 



6t ^ TOFc - TOF^ 

If this time lapse is positive means that the neutrinos arrived faster than expected for light, while if this quantity 
is negative means that the neutrinos are slower than light. 

To obtain the average delay in time, between the neutrino detection and the PDF function, a numerical maximum 
likehood analysis is performed. The results are shown separately for each of the the two extraction, that are enough 
separated by 50ms to be uncorrelated. The results are summarized in figure 1 taken from [l|: 




FIG. 1: Graphic of results of the PDF for protons and the detected neutrinos detected by OPERA 

Surprisingly, a positive lapse shift of 5t — 1048. Sns was obtained, indicating an early arrive than expected at the 
speed of light, there are corrections due to the electronic time tag GPS, UTC, BCT, etc that are summed up to 
(985.8 ± 7A)ns, when all the chain of systemic errors were taken into account. So, there is still an unexplained 
forward lapse shift of: 

(57.8±7.8±^J)ns 

that indicates that neutrinos are 0.25 x 10"** faster than known speed of light in the vacuum. In this result had 
been included the statistical error ±6. 9ns obtained from the maximum likehood analysis, and checked with various 
combination of MonteCarlo simulations. 
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Recently the OPERA collaboration has reported 20 neutrino events obtained for very short proton extraction pulse, 
with a PDF gaussian mean width of only 3ms, in four bunches separated 524ns that is roughly one thousand times 
shorter than former long events. This 20 neutrinos also appears to precede light for 

62.1 ± 3.7ns with a RMS: 16.4ns 

These short extractions has a much lower particle density with only 1.1 x 10^^ protons on target and mesons average 
lorentz factor as large as ^-n = 190. This result of only 20 events has a very high Root Mean Square of RMS := 16.4ns, 
so it will be significant only when taken jointly with the previous 16000 events result, but it has been used to discard 
the shape deformation factor as a possible explanation to the early neutrino arrival. 

III. DISCUSSION. 

In a recent work j?!, it has been proposed an explanation for the neutrinos coming from the "long" 10500ns LHC 
extractions, that apparently arrives to Gran Sasso 57.8ns earlier with respect to the speed of light in vacuum. That 
apparent early arrival could be explained by a shape deforming effect of the neutrino distribution function, with respect 
to the proton distribution function (PDF) due to the stimulated absorption while traveling by the decay tunnel that 
acts as a LASER tube, where other effects as stimulated emission, absorption and spontaneous emission could take 
place. In that work, it was obtained a shift factor for the maximum of the distribution functions is obtained, due to 
the combined stimulated absorption and emission effect: 

Ai - {N2B21 - (1) 

where N2 is the pion or meson density inside the decay tunnel, TViis the muon density, -B2iis Einstein Q stimulated 
emission coefficient, C12 is stimulated absorption coefficient and a the gaussian root mean square of the initial neutrino 
pulse. 

When stimulated absorption dominates over emission, there is negative or backward shift in time for the maxima 
of the distribution functions while traveling through the L = 1km decay tunnel. If only stimulated absorption is 
considered (Beer-Lambert case, B21 = 0) for the long pulse extractions, with A^i IQ^^ /vol (with a tunnel volume 
vol ~ 3000m^) taking a ~ 3300ns (that is also consistent with the tunnel length) with a stimulated absorption 
coefficient: 

3 

C12 ^ 1.65 X 10"*^— = 5.49 X 10"6arn =4> 2.1'dharnl muon 
sec 

will produce the ~ 60ns that had been reported by OPERA. But for the LHC a = 3ns short pulse extractions, with 
a density of 1.1 x IQ^'^ /vol, and the former C12 coefficient will give an insignificant time delay of 

At = 0.0000015ns 

As equation([T|) is strongly dependent on that change for a factor of 10^ (while density changes by an additional 
factor of 60) stimulated absorption or emission could not explain the early arrival of short and long CNGS beam 
simultaneously. 

IV. SPONTANEOUS EMISSION IN THE CONTEXT OF SPECIAL RELATIVITY 

It is not know in which part of the 1000m = 3300ns tunnel the decay of mesons into neutrino plus muon occurs, 
the statistical correlation of thousands of neutrino peak distribution against the PDF function will give the mean 
delay time with high statistical precision when it is assumed that the spatial distribution of the decays in the tunnel 
is random or gaussian, but it could happens that the starting point of neutrinos could be shifted in time or driven 
by stimulated or spontaneous emission-absorption processes, as happens for a pulse LASER traveling through an 
amplifier plasma with an initial population Q. 

The balance equation considering spontaneous emission {A21 coefficient), stimulated emission {B21 coefficient) and 
stimulated absorption (C12 coefficient) proposed in Q: 
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an an _ dNi dNi 
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+ — •^= — + — -^M (4) 



with neutrino density number n. The neutrino velocity is and the the meson and muon velocities are v-^ and 
respectively. All those densities are functions of the tunnel length x and the time t. For short a = 3ns LHC pulse 
extractions the 3300 ns decay tunnel will be almost empty, so quadratic stimulated term could be a million times less 
than spontaneous emission process, so the system of equations reduces to: 



- + -.. = iV2A2, (5) 
dN2 dN2 , 

^ + ^... = ~A^.A., (6) 

The second equation is Lorentz invariant, but the first is no exactly. That is due that for ^ the 4-divergence is 
taken for a 4- volume enclosed by the transverse area of the tunnel times a length equals i^tt-t and multiplied by a time 
interval equals r — (^21)"^ the mean decay life of the meson. Under a Lorentz transformation, the time length scales 
a Lorentz factor while the spatial length scale by the inverse of the Lorentz factor then the volume in which is taken 
^ is invariant. But for equation ([5]), the volume length equals v.t, where the speed of neutrinos is involved instead 
of that of the mesons, also the time interval to be taken is t = (^21)"^, where the mean decay time life is measured 
in the rest or Lab system and not as seen from the meson. To correct this difference in the volumes in which particles 
are destroyed or created a factor v/v^r must be included but at the rest or Lab system this factor is almost one: 
N2 A2i.v/v.,r — N2 A, but equation ([S]) is valid only in the Lab system reference frame. 

So in order to get the correct solution to the process we must go to the meson reference frame where v-^ = Q and 



then the equation system is: 



dn' dn' , m 



(8) 



dt' To 

where Tq = 26ns is the mean life of the Pion in the meson reference frame. 
Equation (jH]) could be easily integrated 

iV^=7V^„(x')exp(--) (9) 

To 

A^^ =7V^oexp(-i^^)exp(-^) (10) 

iCTo To 

where a gaussian shape as in [3] was used and ctq is in the meson reference frame. Equation ([7|) could also be integrated 
in term of the error function, but is easier to boost back to the Lab reference frame. 

> ( i)77r = -^ttTtt ~ -^Itt (11) 

t'^it- v^x)^^ = + ^l) + —]-ii. + 7r- (^2) 

2 Ttt 

X 

Vtt = t-\ , r] = t^x, (13) 

^■n =t- i^t-X (14) 
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where the approximation w 1 has been used at right hand equations. Then the sohition to ([7]) and ([8]) could be 
written in the Lab Reference frame as a function of hght coordinates {i],^) instead of the more exact characteristic 
coordinates (?77r,C7r) 



N2 = N2oexp{^ 



'^)exp(- 



-— )exp( 

2t To 



^2 

n = N20 exp(-— ^)exp( )[1 - exp(- 

2o-^ To 

where a — (To/jtt and_ r = JttTo are in reference to the Lab Frame. 
Proceeding as in 7] looking for the maximum: 



2t' 



(15) 
(16) 



dn 

at 



= 



dn dn ^ 
To 



then, valuating at the tunnel's end [£] 



{t2,L) 



3 X lO^km/ s. Equation could be rewritten as: 



exp(-^)] 



1 / ^ ^ 



(17) 
(18) 



At — t2 — ti and [rj] 



(*2,i) 



(«2,i) 

At + 2L/c where L = 1km and c 



1 

27 



exp[7^](l +exp[— ] 
Zt tc 



(19) 



taking into account that At 
using that (1 + exp[^]) 



-60ns <<2t — 9880ns, then exp[-^] « 1 in order to solve the transcendent equation 
2.95 equation (fT9|) could be solved by 



At 



77rO- 



-65.8ns + 0.0003ns 



(20) 



To 5.97^To 

As the forward time shift, for the 20 short pulse neutrinos obtained for OPERA is —62.1 ± 3.7ns, the theoretical 
obtained At — —65.8ns completely explain it. This striking result make us to ask, if also the 57.8ns early arrival for 
the "long" 10500ns pulse extractions event could be explained by the 200Mhz LHC harmonic, that is superposed to 
the long pulse with a period of T = 5ns, as could be seen in the following picture taken from yy]. 
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FIG. 2: Graphic of results of the PDF for the protons from LHC showing the 200Mz harmonic. 



If we approximate this harmonic in the interval [— f^, -j] by a gaussian function 

1/ 27rt\ , P , 

-(1 + cos— )«exp(-^) 



for t = a we get: a = 1.08ns that is enough to explain a early arrival of only 

At = = —7.3ns 



So apparently, spontaneous emission only, could not explain the early arrival of long CNGS pulses, but as was done 
in Q the remaining 50ns could be fully explained if stimulated absorption (and possible stimulated emission) is taken 
into account. 
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V. CONCLUSIONS 



In conclusion the deformation of the a — Sns short pulses of neutrinos, calculated with relativistic corrections, 
produce an apparent" earlier arrival of the neutrinos, with respect to the speed of light by 



At 



77rO- 



-65.8ns 



explaining the reported earlier arrival of —62.1 ±3. 7ns completely. So the short pulse neutrinos are not superluminal, 
there is a relativistic time shift that produce a deformation of the wave distribution function that shift the maximum, 
so it will arrive earlier, not real particles. Here is a clear prediction: the time shift for short pulses will be linear in 
the Lorentz factor and quadratic in the mean gaussian width a that could easily tested at OPERA with the actual 
infrastructure. 

As a second conclusion, the spontaneous emission process alone could not explain the —57.8ns. But it could be 
fully explained if stimulated absorption-emission processes were taken into account. The calculation performed in Q 
is a very simple and restricted approximation, exact solutions of ([5]) ([3]) and (U) that takes into account relativistic 
correct volumes and real pion, muon and neutrino velocities are required. 

There is a way to probe the existence of these shape deformation processes: to compare the detected neutrino time 
distribution with the muon probability distribution function, instead of the meson probability distribution function, 
because as the muon and the neutrino are created at the same event: they must have the same forward time shift 
so there, must be not difference in time at which intensities maximum are achieved. In [ij the muon PDF was not 
considered, because the simulations show that will give null corrections. That is not true as was shown in this paper. 
Muon PDF may be difficult to analyze because the muon is a very interactive charged particle, but comparison of 
muons PDF with detected neutrinos must be performed in order to establishes if neutrinos travel faster than light or 
stimulated absorption or emission of neutrinos exist. 
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